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We emphasize and demonstrate that, besides using the usual recursion relations involving succes- 
sive layers, generalized Fresnel coefficients of a multilayer can equivalently be calculated using the 
recursion relations involving stacks of layers, as introduced some time ago [M. S. Tomas, Phys. Rev. 
A 51, 2545 (1995)]. Moreover, since the definition of the generalized Fresnel coefficients employed 
does not imply properties of the stacks, these nonstandard recursion relations can be used to cal- 
culate Fresnel coefficients not only for a local but also for a general multilayer consisting of various 
types (local, nonlocal, inhomogeneous, etc.) of layers. Their utility is illustrated by deriving a few 
simple algorithms for calculating the reflectivity of a Bragg mirror and extending the formula for 
the Casimir force in a planar cavity to arbitrary media. 

PACS numbers: 42.25.Bs,42.25.Gy,12.20.Ds,42.50.Lc 



Generalized Fresnel coefficients are basic ingredients 
in the theory of electromagnetic processes and effects in 
layered systems such as light propagation in stratified 
media [l|, Q , molecular fluorescence and energy- transfer 
near interfaces [3[ , (dipole) radiation from multilayers 0- 
spontaneous emission and light scattering at surfaces 
and in planar cavities fiH, the Casimir effect between 
multilayered stacks [lOl - ll 2^ etc. Correspondingly, the 
problem of calculating Fresnel coefficients arises in many 
area of physics: optics, surface physics, spectroscopy of 
multilayers, cavity QED, theory of the Casimir effect etc. 
For systems consisting of few layers, these coefficients 
are standardly calculated using (ordinary) recursion re- 
lations involving successive layers. With increasing num- 
ber of layers, however, this method soon leads to cumber- 
some formulas and thus becomes impractical. Therefore, 
despite the possibility of a polynomial representation of 
the generalized Fresnel coefficients [l3[ , for more complex 
systems these coefficients are conveniently calculated us- 
ing the transfer matrix method 0,0, Hill- 
Based on the definition of the generalized Fresnel co- 
efficients given in Ref. [6j, in our consideration of the 
Green function for a (local) multilayered system Q we 
have used the recursion relations for the reflection (r) 
and transmission (t) coefficients involving stacks of lay- 
ers, which enabled us to write the Green function in a 
simple compact form. For a stack of layers between, say, 
layers j and m (denoted shortly as j /m = j . . .m) these 
recursion relations read [§[ (unless necessary, we omit the 
polarization index q — p,s) 
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where k is an intermediate layer, as depicted in Fig. [TJ 
and where by using the notation rj/ k / m we simply stress 



to which intermediate layer we address. As seen, the 
above recurrence relations look the same as the standard 
ones 0, 0] (to which they reduce in case of a system 
jk/m = jk . . . m), however, this time they generally in- 
volve Fresnel coefficients for stacks of layers. 
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FIG. 1: Stack considered when deriving recursion relations for 
generalized Fresnel coefficients shown schematically. Layers 
j, k and m are local and described by (complex) refraction 
indexes n a (oj) = ^£ a (uj)(ia(u), a = j,k,m, whereas stacks 
between them are unspecified. Arrows indicate propagation 
of a wave incident on the stack. A shifted- z coordinate system 
is adopted, as explained in the text. 

Although Eqs. (JXJ) and ^ have also been known for 
some time in optics of multilayers [3], it seems that the 
possibility of grouping the layers in stacks when calcu- 
lating Fresnel coefficients is not widely recognized. Since 
this method is particularly convenient in some cases, e.g., 
when calculating Fresnel coefficients of periodic media, in 
this work we explicitly derive the above recurrence rela- 
tions starting from the definition of the generalized Fres- 
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nel coefficients [6j (this derivation was omitted in Ref. 
@ ) and demonstrate their equivaience with the standard 
recursion reiations involving successive layers. We use 
this opportunity to emphasize that Eqs. and © 

are also valid for piecewise nonlocal stratified media and 
adapt their derivation accordingly. We illustrate their 
usefulness on a few simple examples including the theory 
of the Casimir effect in (nonlocal) multilayers. 

Consider a linear system consisting of isotropic layers 
0. . .n = 0/n. The electric field of a wave incident on 
the system with the parallel wave vector k can in a local 
layer / be written as 

E ; (k,a;,r) = E ; (k, uj, z) exp(ik • ry), 

E,(k,a;,«) = e+(k)e^£+ + ef{k)e-W"Ef, (3) 



where /3; = ^Jk 2 — k 2 , with ki — niuj/c — Jejjliuj/c. For 
a local system, the coefficients Ef~ are determined by 
matching the field at boundaries of each layer and can 
be entirely expressed in terms of the incident-wave am- 
plitude and the generalized reflection and transmission 
coefficients of the corresponding stacks of layers [1, 0] ■ 

Wc define generalized Fresnel coefficients of a stack be- 
tween two (local) layers as follows |6[ . A reflection coeffi- 
cient r of a stack is the ratio of the reflected to incoming 
wave electric-field amplitude (factors multiplying e's) at 
the corresponding stack's boundary considering only lay- 
ers of the stack, i.e. as if there are no other layers in the 
system. Similarly, a transmission coefficient t of a stack 
is the ratio of the transmitted to incident wave amplitude 
calculated at the corresponding stack's boundaries as if 
there are no other layers present in the system. Note 
that these definitions do not imply any property of the 
stack. Therefore, we can use these coefficients to describe 
wave propagation not only in a local system but also in 
a general system consisting of various (local, nonlocal, 
inhomogeneous, unspecified etc.) layers. In calculating 
these coefficients it is convenient to adopt a (shifted-z) 
representation for the field [9j in which < z < di in any 
finite layer I, whereas — oo < z < and < z < oo in 
external layers I = and I = n, respectively. 

On the basis of the above definition we can straightfor- 
wardly calculate the recurrence relations for generalized 
Fresnel coefficients rj/ m and tj/ m of the stack of layers 
between local layers j and m. Indeed, according to its 
definition, to calculate the reflection coefficient fj/ m we 
consider the wave incident from the layer j upon the sys- 
tem j/m (see Fig. [1]). Then E 1 " = and by definition 
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However, rj/ m can also be calculated by considering 
transmission of the wave to an intermediate (local) layer 
k and its subsequent partial transmission back to the 
layer j. Considering the field at the relevant boundaries 
and using the above definitions of the generalized Fresnel 



coefficients, the amplitudes of the wave are then related 
by the following set of equations 
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Eliminating E k and comparing the ratio EJ / E^ with 
that from Eq. ((4]), we arrive at Eq. ([T]). Similarly, the 
transmission coefficient tj/ m of the stack is by definition 
given by 
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On the other hand, by considering transmission of the 
wave to the intermediate layer k and its subsequent par- 
tial transmission to the layer m, we find that the ampli- 
tudes .E+ and Ef are related by Eq. ([5]), with Eq. (|5aj) 
replaced by 



E + — t, 1 p'^kdk J?+ 
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Proceeding as before, we obtain Eq. ([2]). 

Clearly, Fresnel coefficients must not depend on the 
choice of the intermediate layer in Eqs. ([]} and (|2|). 
Therefore, to prove the consistency of these recurrence 
relations, we must show that 

Tj/k/m ^j/l/m: ^j/k/rn tj/l/rm (p) 

where I denotes some other intermediate local layer. To 
prove this for rj/ m , we rewrite Eq. (JTJ) in the form 
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consider the layer I between layers k and m and apply 
Eq. (|9]) to the reflection coefficient rk/ m — r k/i/m in this 
very same equation. Rearranging the terms and using 
again Eq. © to recognize reflection coefficients rjn and 



r i/j 



we find 
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Now, noting from Eqs. @ and ([S]) that 
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and that fj/ifmJs given by a similar expression, we find 
that actually 



Pi 

El 



a 3/l 



rj/m/j =0,3/1. 



(12) 



Accordingly, the right-hand side of Eq. (TTU)) is, upon 
using Eq. ©, indeed identified as Tj/i/ m . To prove the 
equality of the transmission coefficients in Eq. (|5J) we 
first apply Eq. ^ to the coefficient thjm 
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then use the identity (r k/m = r k/l/ri 

(1 ~ r l/k r l/m e 2 ^ d ')(l - r k/j r k/m e 2i ?" d *) = 
(1 - r k/j r k/l e 2i ^ d ")(l - n/jn/m^*), (13) 

which follows from Eq. ([9]). In this way, we obtain 
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which is according to Eq. ([2]) equal to tj/i/ m . 

In most textbook approaches to the wave propagation 
in layered media the basic ingredients in the calculation 
of the generalized Fresnel coefficients are coefficients rj k 
and tj k for the interface between two neighbouring local 
media j and k. With unit polarization vectors in Eq. ([3]) 



TAk+fcz 



k x z, 
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and applying the usual boundary conditions for the field, 
it is straightforward to show that the above definition of 
Fresnel coefficients leads to the standard single-interface 
coefficients 
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where 7j fc = ej/e k and j? k — jjLj/^i k . As can be easily 
verified, these coefficients obey the Stokes relation 

a jk = tjktkj ~ r jk r kj = 1. (17) 

We note that the symmetry property of the single- 
interface transmission coefficient, as expressed by Eq. 
(|16b[) . implies for local systems the same symmetry prop- 
erty of the generalized transmission coefficient tj/ m , that 
is 11 



f^mPj 



(18) 



Indeed, through Eq. ©, it certainly holds for a three- 
layer system jkm owing to the symmetry property of the 



single- interface transmission coefficients tj k and t km . As- 
suming in Eq. ([2]) the same symmetry properties of the 
generalized coefficients tj/ k and t k / m , we immediately 
find that Eq. (fT8| is obeyed. Accordingly, using the in- 
duction argument, we may conclude that this equation is 
for local systems generally valid. Now we argue that Eq. 
(fl"8]l is also valid for piecewise nonlocal systems as it actu- 
ally ensures the equality of the stack's transmittances for 
waves incident on it from either side. Indeed, the Poynt- 
ing vector of the upward/downward-propagating wave in 
a layer I is according to Eqs. Q and (fT5"|) given by 



P±(k,u,,z) = ±ReVm^^\Efe^\ 2 , (19) 
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with = Tjf* = £*/fJ-i- Accordingly, assuming the out- 
most layers j and m transparent, transmittances of the 
j/m stack for waves incident upward and downward on 
it [given by the ratios of the respective transmitted- to 
incident-energy fluxes 7}/ m = P+ Z (k, 0)/P+ (k, w, 0) 
and T m /j = Pjz (k, u, 0)/P mz (k, oj, 0)] 



<r _ ^jPm i , i 
I j/m — a \ b j/m\ 



%n/j — 
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can only be equal if Eq. ([18]) is fulfilled. Thus, being 
a consequence of the reciprocity property of the electro- 
magnetic field, this equation is valid for a quite general 
class of (nongyrotropic) media [l6j |. 

As follows from the above results, recurrence relations 
for Fresnels coefficients of a multilayered system can be 
written in a number of ways depending on number of in- 
termediate local layers. We illustrate this by calculating 
Fresnel coefficients of few simple systems. We first con- 
sider a 12/3 system (that is, we allow for an unspecified 
stack of layers between layers 2 and 3 of a three-layer 
local system 123). Since from Eq. (TTT]) a 12 = 1, Eq. © 
leads to standard forms of the recurrence relations for 
Fresnel coefficients 21 
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ri2 + r 2 / 3 e 



2i/3 2 d 2 
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tuti/se^ 
- r2ir 2 /3e 2l/3 



(21) 

which, for a three-layer system 123, reduce to the well- 
known results usually quoted in textbooks [l|, Q. Next, 
we consider Fresnel coefficients of a 123/4 system. Ac- 
cording to Eq. ([9]), they can be calculated from two 
equivalent sets of recurrence relations differing in the 
choice of the intermediate layer. The first set is the stan- 
dard one Q and is given by Eq. ([2"Tj) . with the replace- 
ment of indices 2/3 —> 23/4. The second set of recurrence 
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relations for r 12 3/4 and £123/4 reads 
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where 7-123, r 32 i, £123 and t 32 i are given by Eq. (|2"Tj). In 
our consideration of the Casimir effect in a planar cavity 

llj we have used this result to calculate the reflection 
coefficient of a slab in front of a cavity mirror. 

As mentioned, grouping the layers is particularly useful 
when calculating Fresnel coefficients of periodic media. 
Consider, for example, the reflection coefficient of the 
central segment 121 . . . 121 of a Bragg mirror formed by 
two alternating quarter- wavelength (d, = A/4rij) dielec- 
tric layers 1 and 2. Noting that at normal (k = 0) inci- 
dence e 2ll3id ' — —1 and am = —1, we find from Eq. (|22|) 
(with 3 and 4 — > 1) that the normal reflection coefficient 
Rn of the segment with TV type 2 layers can be calculated 
using the algorithm: Rn = {Ri + i?7v-i)/(l + R\Rn-i), 
with i? = and i?i = r V2 \ = 2ri 2 /(l + r\ 2 ) as in- 
put. Of course, grouping the layers in a different way 
leads to a different algorithm. Thus, one may show 
that starting with R 2 the following algorithm holds: 
R4N — 2R 2 n/(1 + R 2N ), where in each step the num- 
ber of type 2 layers is doubled [17 1. 

Evidently, the above nonstandard recurrence relations 
are particularly convenient when the Fresnel coefficients 
of a stack are already known (being either calculated sep- 
arately or measured) or are to be calculated by some 
method at a later stage. To illustrate this, we (re)derive 
and generalize to media with arbitrary properties the for- 
mula for the Casimir force on a slab in a planar cavity 

11]. Referring to Fig. [5] for the description of the sys- 



tem, the Casimir force (per unit area) acting on the slab 
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FIG. 2: A slab (s) in a planar cavity (ni = 712 = n c ) schemat- 
ically. Cavity mirrors are described by their reflection co- 
efficients Ri and R2 and the slab by its Fresnel coefficients 



ri/2 = r 2 /i and t = t 
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with Tz£ being the relevant component of the vacuum- 
field (Minkowski) stress tensor in the cavity region j. 
Here n = \J n 2 (i^)^ 2 /c 2 + k 2 is the perpendicular wave 
vector at the imaginary frequency (uj = i£) in the cavity 
and rj±(i£ > , k) are the reflection coefficients of the right 
and left stack of layers bounding the region j. We ob- 
serve that according to the identity Eq. (fT3|) the tensor 
components Tj$ and are related to each other [ill ]. 

Noting that ^i-( 2 +) = -^1(2) an( l using the recurrence 
relation [cf. Eq. {3}] 



n+(2-) = r + 



t 2 R 2{1) e- 2Kd ^) 
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Since properties of the slab and mirrors are not specified 
this result is valid for arbitrary media. For local uniform 
media described by the corresponding refraction indexes, 
Fresnel coefficients r — r\ s2 and t — t\ s2 are given by Eq. 
(f2"Tj) and those of the mirrors by Eq. p^|) . In that case, 
of course, the above result agrees with the one obtained 
through a conventional way [l8j ■ When removing a mir- 
ror (letting, say, d 2 — > 00), the above result gives the 
Casimir force between two arbitrary planar objects, as 
also obtained recently through a different approach [19j . 

To summarize, we have emphasized and demonstrated 
the possibility of grouping the layers into stacks when 
calculating Fresnel coefficients of a multilayered system. 
This enables one to consider wave propagation in local as 
well as in piecewise nonlocal stratified media on an equal 
footing. As an example, we have shown that the formula 
for the Casimir on a slab in a planar cavity derived con- 
sidering local media [ll[ is also valid when the objects 
involved have arbitrary properties. 
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